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Abstract 


Wave  propagation  and  random  media  are  defined  and  the  nature  of  the 
mathematical  problems  arising  in  wave  propagation  in  random  media  is  described. 
The  two  principal  types  of  methods  for  solving  these  problems  -  honest  and 
dishonest  methods  -  are  explained.  These  methods  are  first  illustrated  by 
considering  the  geometrical  optics  of  a  random  medium  by  one  method  of  each 
type.  Some  nev  results  are  obtained  by  an  honest  method  and  some  errors  in 
a  previous  work  are  pointed  out.  Comparison  is  made  between  the  results  of 
the  two  methods  and  the  reasons  why  they  disagree  are  explained.  As  a  sec¬ 
ond  illustration  of  an  honest  method,  an  analysis  of  the  reduced  wave  equa¬ 
tion  in  a  random  medium  is  presented.  Some  known  results  are  obtained  in  a 
new  way  which  is  simpler  than  the  usual  one  and  which  appears  to  be  capable 
of  yielding  further  results. 


Table  of  Contents 

Ises. 


1.  Introduction  1 

2.  An  application  5 

3«  Light  rays  in  a  slightly  inhomogeneous  medium  7 

4.  light  rays  in  a  random  medium  9 

3.  Light  rays  as  a  Markoff  process  13 

6.  Phase  fluctuations  20 

7.  Amplitude  fluctuations  2k 

8.  The  reduced  wave  equation  in  a  random  medium  27 

35 


References 


1. 


Introduction 


Wave  propagation  is  one  of  the  means  by  which  energy  travels.  'Propa¬ 
gation'  is  the  process  whereby  the  energy  moves  from  one  region  of  space  to 
an  adjacent  region,  and  'wave'  is  a  general  term  for  a  moving  spatial  dis¬ 
tribution  of  energy.  The  matter  in  the  region  of  space  through  which  the 
propagation  occurs  is  called  the  transmission  medium.  A  random  or  stochas¬ 
tic  medium  Is  a  family  of  media  together  with  a  probability  distribution 
which  specifies  the  probabilities  of  the  various  members  of  the  family. 

Thus  wave  propagation  in  a  random  medium  refers  to  propagation  in  each  mem¬ 
ber  of  the  family  of  media,  together  with  the  probability  of  each  member. 

This  probability,  when  associated  with  the  wave  motion  in  each  medium, 
characterizes  a  random  wave  motion. 

Mathematically  a  wave  motion  is  described  by  a  vector-valued  function 
u(x,t)  of  the  position  vector  x  and  the  time  t.  As  a  consequence  of  the 
physical  lawB  governing  the  wave  motion,  the  function  u(x,t)  satisfies  cer¬ 
tain  equations.  Usually  they  are  partial  differential  equations  pf  hyper¬ 
bolic  type  and  often  of  symmetric  hyperbolic  type.  The  transmission  medium 
is  characterized  by  a  vector- valued  function  n(x,t)  which  enter  the  coeffi¬ 
cients  of  the  equations.  A  wave  propagation  problem  is  that  of  determining 
a  solution  u(x,t)  of  the  equations  which  satisfies  certain  auxiliary  condi¬ 
tions.  These  conditions  are  usually  initial  and  boundary  conditions.  The 
problem  is  Baid  to  be  well  set,  well  posed  or  properly  posed  if  it  possesses 
a  unique  solution  which  depends  continuously,  in  an  appropriate  norm,  upon 
the  data  given  in  the  auxiliary  conditions  and  upon  the  coefficients  in  the 
equation. 

A  random  medium  is  characterized  by  a  family  of  functions  n(x,t,u>)  de¬ 
pending  upon  a  parameter  w.  Here  u  is  a  point  In  a  probability  space  Q  in 
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which  a  probability  measure  dP(to)  is  defined.  Of  course  the  function 
n(x,t,w)  must  be  measurable  with  respect  to  dP(ui).  If  a  well-posed  wave 
propagation  problem  is  formulated  for  the  random  medium,  it  has  a  solu¬ 
tion  u(x,t,w)  for  each  to.  This  family  of  solutions  describes  a  random  wave 

motion.  The  equations  satisfied  by  u  with  the  coefficient  n(x,t,w)  are 

« 

called  stochastic  equations  -  stochastic  differential  equations  if  they  are 
differential  equations.  Thus  in  this  case  wave  propagation  in  a  randum  me¬ 
dium  is  part  of  the  subject  of  stochastic  differential  equations. 

The  reasons  for  studying  random  wave  propagation  are  exactly  the  same 
as  those  for  studying  any  other  random  phenomenon.  There  are  essentially 
two  such  reasons.  In  the  first  place,  we  may  wish  to  consider  a  case  in 
which  the  medium  is  not  known  precisely,  but  in  which  it  is  known  to  be  a 
member  of  a  certain  family  of  media.  If  we  also  know  the  probability  that 
it  is  any  particular  member  of  the  family,  we  can  determine  the  probability 
that  the  wave  motion  is  any  one  of  the  associated  family  of  wave  motions. 
Then  we  can  determine  the  expected  or  mean  wave  motion,  as  well  as  its 
variance  and  other  statistics.  This  statistical  information  can  be  used 
to  estimate  what  is  likely  to  be  the  wave  motion  in  the  case  under  consider¬ 
ation.  In  particular  some  statistical  properties  of  u  may  depend  only  upon 
known  statistics  of  n.  When  we  consider  many  media,  all  of  which  are  mem¬ 
bers  of  the  same  family,  we  may  expect  the  statistics  of  the  corresponding 
wave  motions  to  correspond  to  those  of  the  random  wave  motion  occurring  in 
the  random  medium. 

In  the  second  place  we  may  wish  to  consider  a  known  medium  with  very 
complex  properties.  Then  the  associated  wave  motion  will  also  be  very  com¬ 
plex  and  its  precise  deteimination  may  be  impractical.  Certain  statistical 
properties  of  the  wave  motion  in  a  random  medium,  of  which  the  medium  in 
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question  is  a  member,  may  be  less  complex  and  more  easily  determinable.  In 
other  words,  it  may  be  easier  to  solve  the  random  problem.  If  the  random 
medium  is  appropriately  chosen,  these  statistical  properties  may  be  closely 
related  to  the  actual  properties  of  the  vave  motion  in  the  known  medium 
under  consideration. 

We  shall  see  that  there  are  two  kinds  of  methods  used  to  solve  problems 
of  wave  propagation  in  random  media.  These  may  be  called  'honest'  and  'dis¬ 
honest'  methods  respectively.  In  an  honest  method  the  solution  u(x,t,u)  is 
first  determined  for  each  value  of  u>.  !Ms  solution  may  sometimes  be  found 
exactly  and  explicitly,  but  more  often  it  is  expressed  in  the  form  of  a 
series  in  some  parameter,  or  els  a  sequence  of  iterates,  or  by  some  other 
approximation  procedure.  In  the  process  of  solving  for  u(x,t,u>)  randomness 
plays  no  role  and  therefore  it  provides  no  advantage.  The  second  step  in 
an  honest  method  is  to  compute  the  mean  value  of  u(x,t,u>),  as  well  as  its 
variance  and  other  statistics,  from  the  explicit  expression.  In  this  step 
randomness  may  have  the  helpful  effect  of  yielding  simpler  expressions  for 
the  statistics  of  u  than  those  for  u  Itself. 

In  a  dishonest  method  randomness  is  utilized  before  u(x,t,u)  is  deter¬ 
mined.  For  example,  if  the  mean  value  <u>  is  sought,  the  original  equations 

for  u  may  be  'averaged'  to  yield  equations  for  <U>.  In  such  cases  the  aver- 

2 

aged  equations  also  involve  terms  like  <U  >  or  <nu>.  Dishonesty  enters 

2  2 

through  the  assumption  that  <U  >  ■  <u>  or  <3iu>  =  <HXU>  or  some  similar 
unproved  assumption.  The  reason  for  making  such  assumptions  is  to  obtain 
determinate  equations  for  <u>.  Other  examples  of  dishonest  methods  will  be 
described  later.  In  all  cases  probability  is  introduced  before  u  is  deter¬ 
mined  and  an  improved  assumption  is  made  about  some  statistical  property  of 
the  random  wave  motion.  The  assumption  usually  simplifies  the  problem  so 
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that  it  becomes  solvable. 

Dishonest  methods  have  the  advantage  over  honest  ones  that  they  simplify 
the  problem  to  be  solved.  As  a  consequence,  a  problem  which  can  be  solved 
honestly  only  by  a  perturbation  method  might  be  solved  dishonestly  without 
the  use  of  a  perturbation  expansion.  Then  the  dishonest  solution  would  be 
applicable  for  all  values  of  the  relevant  parameter  while  the  honest  solution 
would  be  valid  only  for  small  values  of  it.  For  this  and  similar  reasons, 
many  of  the  significant  and  non-trivial  results  in  the  theory  of  wave  propa¬ 
gation  in  random  media  have  been  obtained  by  dishonest  methods.  Many  of 
these  results  have  compared  well  with  experiment.  Thus  one  of  the  important 
mathematical  problems  in  this  field  is  to  Justify  the  dishonest  methods  by 
showing  that  their  results,  in  some  sense,  are  approximations  to  honest  solu¬ 
tions.  A  clear  understanding  of  the  circumstances  in  which  this  is  the  case 
would  permit  the  introduction  of  many  more  useful  dishonest  procedures. 

Most  of  the  work  on  this  subject  has  been  done  since  19^5/  having  been 
stimulated  by  practical  problems  of  radio  wave  propagation  through  the  atmos¬ 
phere  and  ionosphere,  sound  wave  propagation  in  the  ocean  and  the  atmosphere, 
light  transmission  through  the  atmosphere,  etc.  The  recent  book  by  L-A. 
Chernov M  contains  a  rather  complete  bibliography  and  an  understandable 


account  of  the  present  state  of  this  subject.  Additional  material  is  con- 

M 

talned  in  the  book  by  V.I.  Tatarski L J.  The  related  subject  of  random  wave 
motion  in  a  non-randan  medium  is  surveyed  by  M.  Born  and  E.  Wolf  Chap’  1C^ 

[V 

Some  other  work  is  contained  in  the  symposium  volumes  edited  by  W.C.  Hofftaan1-- 
and  Z.  Kopal  H  .  Because  of  the  complete  bibliographies  in  these  books,  we 


shall  give  relatively  few  references. 
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Physicists  have  studied  extensively  the  propagation  of  vaves  through 
random  collections  of  discrete  scatterers.  A  clear  formulation  of  a  prob¬ 
lem  of  this  type,  together  with  a  new  dishonest  method  for  treating  it,  was 
given  by  L.L.  Foldy ^ .  This  method  was  generalized  to  other  problems  by 
M.  LaxUJ,  who  has  also  reviewed  much  of  the  previous  work  in  this  field. 

Many  problems  have  been  treated  by  this  method  and  others  by  V.  Twersky®  . 
Recently,  J.  Bazer^  proved  that  for  a  one -dimensional  scattering  problem, 
this  method  does  yield  the  correct  result.  I.  Kay  and  R.A.  Silverman 
investigated  the  accuracy  of  another  method,  the  single  scattering  approxi¬ 
mation,  by  determining  the  extent  to  which  randomness  reduces  the  importance 
of  multiple  scattering. 

The  next  section  describes  how  wave  propagation  in  a  random  medium 
applies  to  the  problem  of  the  twinkling  of  a  star.  In  Section  3  we  consi-. 
der  light  rays  in  a  slightly  inhomogeneous  medium,  and  in  Section  b  we  apply 
these  results  to  a  random  medium  obtaining  some  new  results.  We  then  compare 
these  honest  results  with  corresponding  dishonest  ones  in  Section  5,  where 
we  treat  a  light  ray  in  a  random  medium  as  a  Markoff  process.  In  Sections 
6  and  7  we  use  the  results  of  Sections  3  and  4  to  determine  the  phase  and 
amplitude  fluctuations  of  a  wave  in  a  random  medium.  Our  results  are  exten¬ 
sions  of  previously  known  oneB.  This  completes  our  discussion  of  geometri¬ 
cal  optics  in  a  random  medium  .  In  the  final  section  we  present  a  brief  new 
treatment  of  the  reduced  wave  equation  in  a  random  medium.  This  section  is 
independent  of  the  preceding  ones. 


2.  An  application 

Before  illustrating  the  techniques  vised  to  analyze  wave  propagation  In 

f 

random  media,  we  shall  describe  a  physical  phenomenon  in  which  such  propaga¬ 
tion  plays  a  role.  It  is  the  scintillation  or  twinkling  of  a  star.  The  most 
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appropriate  theory  of  wave  propagation  to  describe  this  phenomenon  is  geo¬ 
metrical  optics.  According  to  this  theory,  light  travels  along  certain 
straight  or  curved  paths  called  rays.  These  rays  are  determined  by  ordinary 
differential  equations  in  which  appears  the  index  of  refraction  n(x,t),  a 
scalar  function  which  characterizes  the  transmission  medium.  This  medium  is 
the  earth's  atmosphere  in  the  present  case.  If  n=I  the  rays  are  straight  lines 
which  emanate  from  the  star  in  all  directions.  One  of  them  (really  a  narrow 
beam)  enters  the  eye  of  an  observer  who  is  viewing  the  star.  The  direction 
from  which  the  ray  enters  the  eye  is  the  apparent  direction  of  the  star. 
Actually  n(x,t)  differs  from  unity  by  a  small  amount.  As  a  consequence,  the 
rays  deviate  slightly  from  straight  lines  and  enter  the  observer's  eye  from 
slightly  different  directions  at  successive  instants  of  time.  Therefore  the 
star  appears  to  be  moving  about  its  mean  position.  Its  apparent  intensity 
also  fluctuates. 

We  would  like  to  calculate  the  apparent  direction  and  intensity  as  func¬ 
tions  of  time.  To  do  so  it  would  be  necessary  to  know  the  index  n(x,t).  This 
is  practically  impossible  because  the  variation  of  n( x, t )  with  time  and  posi¬ 
tion  results  from  the  trubulent  motion  of  the  atmosphere.  Consequently  analy¬ 
sis  would  appear  to  be  impossible.  In  the  face  of  this  difficulty  we  treat 
the  atmosphere  as  a  random  medium.  The  random  medium  must  be  so  chosen  that 
its  important  statistical  properties  correspond  to  measurable  properties  of 
the  atmosphere.  In  this  way  we  give  up  the  possibility  of  calculating  the 
apparent  direction  and  intensity  at  any  particular  time.  Instead  we  can  calcu¬ 
late  statistical  properties  of  the  apparent  direction  and  intensity  which  may 
be  related  to  the  actual  temporal  distribution  of  apparent  directions  and  in¬ 
tensities.  We  shall  consider  this  example  further  in  the  next  section.  In 
doing  so  we  shall  make  use  of  the  fact  that  a  ray  traverses  the  atmosphere  so 
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quickly  that  the  Index  does  not  change  significantly  during  the  traversal. 
Consequently  ve  may  assume  that  the  index  is  independent  of  time  in  deter¬ 
mining  the  rays. 


5.  Light  rays  in  a  slightly  inhomogeneous  medium 

Let  n(x, e)  denote  the  index  of  refraction  of  a  random  medium,  which  may 
he  written  as 

n(x,s)  =  1  +  e  n(x)  .  (1) 


The  quantity  €  measures  the  deviation  of  the  index  from  unity.  The  index  also 
depends  upon  a  parameter  to  which  ve  shall  not  write  explicitly.  We  wish  to 
determine  the  ray  x(s,e)  which  starts  from  the  origin  in  the  direction  of  the 
unit  vector  u.  Here  s  denotes  arclength  along  the  ray.  The  ray  also  depends 
upon  w.  The  equations  which  x  satisfies  are 


(nx1 )  «  Vn 

(2) 

XI 

0 

N 

O 

(5) 

x'(0)  *»  u  ,  (u2  »  1)  . 

(*) 

These  equations  have  a  unique  solution  which  depends  continuously  upon  n  and  u 
in  an  appropriate  norm,  if  n  is  continuously  differentiable.  Thus  the  problem 
(2)  -  (4)  is  well  posed. 

To  find  x  we  shall  determine  its  derivatives  with  respect  to  e  at  £  c  0 
and  then  express  x  by  means  of  its  Taylor  series  in  e.  Thus  we  shall  employ 
an  honest  method,  in  the  first  phase  of  which  probability  plays  no  role.  Let 
us  first  set  e  »  0  in  (2)  -  (4)  and  denote  x(s,0)  by  xq(s).  Then  we  obtain 


x"  =  0 
o 

xq(0)  -  0 


(5) 

(6) 

(7) 


m  U 
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The  solution  of  (5)  -  (7)  is 


x0(s)  -  ^  .  (8) 

Nov  ve  differentiate  (2)  -  (4)  with  respect  to  e,  set  e  =  0,  and  denote 

x  (s,0)  by  x  (s),  obtaining 
€  € 


“  ^  n(x J  - 


V  n(xQ) 


(9) 


x#(0)  =  X£(0)  -  0  .  (10) 

— ►  — ► 

The  right  side  of  (9)  is  Just  the  component  of  Vn  which  is  normal  to  x'  =  u. 
Let  us  call  it  the  transverse  gradient  and  denote  it  by  .  Then  the  solu¬ 
tion  of  (9)  and  (10)  is 
s 

x€(s)  »  J  (s-t)  VT  n(ut)dt.  (11) 

o 

Differentiating  (2)  -  (4)  twice  with  respect  to  e  at  e  ■  0  yields 


xM(s)  =  2(x  •  V  )Vm  (i(x  )-Vn1  fx  (x  )-2x'x'  •  V  n(x  )-2x'x'  •  V  |x(x  )  (12) 

€€ '  €  T  T  '  o  T  o  o  e  T  o  e  o  '  o 


x  (0)  »  x'  (0)  -  0  . 
€€ '  '  €€ 


(13) 


The  solution  of  (12)  and  (13)  is 


x£e(s)  ■  2  J  (s-t)  a 


Xe  ’  VT 


"  I  VT  ^“2(u  •  Vn)Vr 


T 


dt. 


w 


o 

—¥ 

In  the  integrand,  the  argument  of  n  is  ut.  Thus  to  the  second  order  in  e  we 
have 


x(s,e)  =  xQ(s)  +  e  x6(s)  +  x€£(s)  +  0(e5)  . 


(15) 
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4.  Light  rays  In  a  random  medium 

Let  us  now  take  account  of  the  randomness  and  compute  some  statistical 

— ►  , 

properties  of  x(s,e).  Let  us  begin  with  the  mean  value  <k(s,e)>  which  is 
just  the  sum  of  the  mean  values  of  the  terms  on  the  right  side  of  (15). 

Since  x  (s)  =  us  is  independent  of  p.  and  therefore  of  to ,  <xq(s)>  =  us.  From 
(ll),  by  interchanging  the  order  of  taking  the  mean  value  with  integration 
and  differentiation,  which  we  assume  to  be  permissible,  we  obtain 
s 

<Xg(s)>a.  J' ( s-t)  <p(ut)>  dt.  (16) 

o 

In  (l6)  only  <L>  occurs.  We  see  that  it  is  not  yet  necessary  to  specify  the 
probability  space  fi,  the  measure  dP(io)  nor  ji(x,u)  if  instead  we  give  <p(x)>. 
We  therefore  assume  that 

<+i(x)>  a  0  .  (IT) 

From  (16)  and  (17)  it  follows  that 

<x£ ( s ) >  «  0  .  ( 18 ) 

To  compute  the  mean  value  of  x^(  s ) ,  which  is  given  by  (l4),  we  must 
know  the  mean  values  of  certain  quadratic  expressions  in  |j.  and  its  deriva¬ 
tives.  They  can  be  determined  from  the  correlation  function  <4.,.(x1)|i(xg)> 
which  will  also  be  required  in  later  calculations.  Therefore  we  assume  that 
it  is  given  by 

<„(X1)P(X2)>  »  <L2>  N(  I^-Xj  )  .  (19) 

2 

Here  the  mean  square  fluctuation  <h  >  is  assumed  to  be  a  constant  and  the 
correlation  coefficient  N  is  assumed  to  be  a  function  of  the  distance  Jx^Xgl 
only.  By  the  Schwartz  inequality  it  follows  that  |n|  £  1.  Of  course  N(0)=1. 


10 


We  also  require  that  N'(0)  *  0  and  that  N(»)  =  0.  It  follows  that  NrT,(0)  ^  0. 
In  addition  N  should  become  small  beyond  a  distance  a  which  we  shall  call  the 
correlation  distance. 

2  2  2 
From  the  assumption  that  <|i  >  is  constant  we  have  <7^  |i  >  =  V^,<h  >  =■  0. 

Thus  the  mean  of  the  second  term  in  the  integrand  in  (l4)  is  zero.  To  com¬ 
pute  the  mean  of  the  first  term  we  use  (ll)  for  and  the  first  term  becomes 

s  *2  r 

/ /  (s-t2)(t2-t1)  J  [vT1^ut1)  •  VT2J  VT2n(ut2) 

0  0  l 

The  expression  in  brackets  can  be  rewritten  as  (V^  *  VT2^VT2M^^t2^^tl)  811(1 

since  V  „  -  V  ,  when  applied  to  functions  of  t0-  t  its  mean  value  Is 
T1  T2  ^  -1- 

<V?>(  VT1  •  VT2) VT2N(|t2-t1j )  -  -  <4i2>VT2  [v|2N(  I  t2-t1 1 )  •  (21) 

2 

If  we  set  | t2-t1|  =»  r  then  N  -  N(r)  and  VT2  =  VT.  Since  VT  N(r)  is  a  func¬ 
tion  of  r  only,  its  transverse  gradient  Is  zero.  Therefore  (21)  and  the  mean 
of  (20)  are  zero.  In  a  similar  way  we  find  that  the  mean  of  the  last  term  in 
(l4)  is  also  zero.  Thus  only  the  third  term  remains  and  we  obtain 

s 

<5c(s,c)>“  us  -ue2  J  (s-t)<(VT  |i)2>  dt  +  O(s^)  (22) 

o 

To  evaluate  <(VT  |i)2>  we  have,  since  VT1  N  =  -  VT2  N  and  since  V^N(O)  = 

|  V^(0), 

<(V  ^)2>  -  <U2>  lim  V  •  VT2  N(  1%-Xgl) 
xx-  x2 

-  -  <U2>  v2n(o)  -  -  |  <H2>  v2n(0)  .  (25) 

Thus  finally*  taking  account  of  the  fact  that  V  N(0)  «  5Nrr(0),  w  e  obtain 
<x(s,c)>  -  u(s  +  €2  <H2>  nrr(0)s2)  +  0(€5)  .  (24) 
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From  ( 24-)  we  can  solve  for  the  length  s  such  that  a  ray  of  this  length 
travels  a  mean  distance  L  in  its  original  direction.  This  value  of  s  is 


s  «  L  -  e2  <ji2>  Nrr(0)L2  +  O(s^)  .  (^5) 

Equation  (2k-)  gives  the  mean  location  of  the  endpoint  of  a  ray  of  length 
s  which  starts  from  the  origin  in  the  direction  u.  This  mean  position  is  on 
the  straight  line  through  the  origin  in  this  direction,  which  we  expect  to  be 
the  case  from  symmetry.  However,  its  distance  from  the  origin  is  less  than  s, 
as  we  also  expect.  The  terms  given  in  (24)  do  not  suffice  to  compute  <Sc(s,€)> 
for  very  large  values  of  es.  This  can  be  seen  by  noting  that  for  such  values, 
the  mean  position  of  the  endpoint  would  lie  in  the  direction  -u,  according  to 
(24)  since  N  (0)  S  0.  This  result  is  unreasonable,  as  is  the  fact  that  the 
mean  distance  from  the  origin  to  the  endpoint  would  exceed  s  for  very  large 
€S.  We  conclude  that  the  approximation  (24)  for  <5c(s,€)>  is  not  uniform  in  s. 
Thus  more  terms  in  the  series  are  required  for  large  values  of  es.  They  in¬ 


volve  third  and  higher  order  correlation  functions  and  moments  of  |i. 

Let  us  now  compute  the  mean  value  of  [X(S,€)]  .  From  (15)  we  have 

J^(s,e)  *  b2  +  e2  x€(s)  +  e2  su*x£e(s)  +  O(e^)  .  (26) 

In  writing  (26)  we  made  use  of  the  fact  that  u  •  X£  *  0.  By  using  (ll)  and 
(l4)  in  (26)  we  obtain 


The  mean  value  of  the  third  term  on  the  right  side  in  (27)  has  been  computed 
in  reducing  (23)  to  (24).  To  compute  the  mean  of  the  second  term  we  write  it 


as  follows 
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»  b 

<(xg)2>=<  /(s-VV^dt,./  (s-t2)VT(ut2)dt., 


s  s 


<^2>  y J  (s-t1)(s“t2)VTi  *  VT2  N(  l^i"  "tglJdt^tj 


0  0 


s  s 


-  -  <M2>  //(s-t1)(s-t2)  vT2  HCjt^tgDat^tg  . 


o  o 


We  now  define  r  and  r  by 

o  J 


tl  -  r0  +  r/2  *2  “  r0  "  r/2  • 

Then  (28)  becomes,  if  we  define  N(-r)  «  N(r), 


o  fl+(r/2) 

r 

o] 

f  1 

t  \2  r 

(s-ro)  -  T 

s  -p/2 

— j 

S  8- 

■(r/2) 

- 

-  <h2>  f 

f 

(s-r0) 

0 

r/2 

VT  W(r)drQdr 


,2  r 


VT  N(r)drQdr 


Upon  performing  the  integration  with  respect  to  rQ  and  combining  the 


integrals  we  obtain  from  (30) 

s 

<(x€)2>»  -  <n2>  J  V2N(r) 

o 

To  simplify  (jl)  we  note  that 


r5  2  2s5 

T-rs  +  -y- 


dr 


V2  N(r)  =  2r-1  »r(r) 


(28) 

(29) 

(30) 

resulting 

(31) 


(32) 
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With  the  aid  of  (32),  ve  can  rewrite  (31)  as 


:2  <#/>„  -  e2  <u2>  J' Ny(r) 


2  r 
3 


dr 


(33) 


Now  (27)  becomes 


<|x(s,e)  |2>  *.  s2  -  e2  <p2> 


r-  B 


/  V 


r)<  “  28 2+  ^  “  2s5  Nrr(0) 


0 


+0(£5) 


w 


For  values  of  s  large  compared  to  the  correlation  distance  a,  (3^)  can  be  sim¬ 
plified  to 


<|x(s,c)  j2>  «  s2-  e2  <n2> 


\  f  r_1  Nr(r)dr  -  2Nrr(0) 


o 


+  0(e5)..(35) 


Equation  (34-)  or  (35)  gives  the  mean  square  distance  between  the  end¬ 
points  of  a  ray  of  length  s.  Let  P  denote  the  transverse  displacement  of  the 
endpoint  of  a  ray  of  length  s  from  the  straight  line  tangent  to  the  ray  at  its 
initial  point.  Then  equation  (33)  gives,  to  order  c  ,  the  mean  square  value 


<p  >.  For  a  large  compared  to  the  correlation  distance  a,  (33)  yields  for 

^€2  S5  ” 


<P2> 


2 

<p  >  m  - 


J  r  1  Nrdr  +  O(e^)  . 


(36) 


By  substruction  of  (36)  from  (35)  we  obtain  the  mean  square  distance  which 
the  ray  travels  along  its  original  direction.  For  s  large  compared  to  a,  it 
is 

<|u*x(s,€) |2>  »  s2  +  2e2  <p2>  Nrr(0)s^  +  O(e^)  .  (37) 
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Let  us  finally  calculate  the  mean  square  value  of  x'(s,e)  -  u.  This  is 
the  mean  square  value  of  2  -  2  cos  a(s,e)  where  a  is  the  angle  between  x'  and 
u.  From  (15)  we  have 


/-*•.  -^2  2/-».k2  •  f  3\ 

(x*  -  u)  -  e  (xp  +  0(0  . 

From  (38)  and  (11)  we  obtain 

<(  x'  -  u)2>  a  e2  <  (^J  VT  n  dt^  >  +  0(e5)  . 

By  proceeding  as  in  (28)  -  (33)  we  find 


(38) 


(39) 


<(x*  -  u)2>*  2  -  <cosa(s,£)>^  »  4fi2  <ja2~ 


N(s)-l-s  r  dr 


+  0(e5). 

(40) 


For  s  large  compared  to  a  this  becomes 


00 

<(x'  -  u)2>-  2^1  -  <cosa(s,€)>)  *5  -  4g2  <p2>  s  J  r-1Nrdr  +  0(e5)  . 


(41) 


If  the  right  side  of  (4l)  is  small  compared  to  unity,  we  have 

2 

2(1  -  <cos  a>)  v  <a  >.  Then  (4l)  gives  the  mean  square  value  of  a. 

The  results  (36)  and  (4l)  can  be  written  simply  in  terms  of  the  ray  dif¬ 
fusion  coefficient  D  which  is  defined  by 


D  ■  -  € 


<p2>  J  r  ^  Nr(r)dr  . 


Then  (36)  and  (4l)  become 

<p2>  «  ^  Ds5  +  0(e5). 

2  5 

<x  >  »  4ds  +  0(c  ) . 


(42) 


(43) 

(44) 
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An  incorrect  form  of  (4l)  vas  derived  by  L.A.  Chernov  [l,p.l7,eq. Jl] • 

His  method  of  derivation  is  the  same  as  that  given  here  but  he  did  not  intro¬ 
duce  a  small  parameter  e  and  expand  systematically  in  powers  of  it.  As  a 
consequence  he  made  an  error  in  going  from  his  eq.  (25)  to  eq.  (26)  on  page 
16.  Therefore  his  result  eq.  (51)  and  his  diffusion  coefficient  eq.  (3i0 
page  17  are  incorrect. 


5.  Light  rays  as  a.  Markoff  process 

In  the  preceding  section  ve  treated  light  rays  in  a  slightly  inhomogene¬ 
ous  random  medium  by  an  honest  method.  We  computed  the  mean  values  of  some 
quantities  associated  with  the  rays  and  obtained  simple  expressions  for  them 
when  the  ray  length  s  was  large  compared  to  the  correlation  distance  a.  Al¬ 
though  the  results  are  correct,  they  are  not  useful  for  very  large  values  of 
s,  as  can  be  seen  by  examining  them.  We  must  conclude  that  the  quantities  we 
have  calculated  in  powers  of  e  are  not  uniform  in  s. 

To  obtain  results  which  are  valid  for  large  values  of  s  a  dishonest  method 

Til  - 

has  been  employed1--1.  It  consists  in  treating  the  ray  direction  x'(s)  as  a 
random  function  of  the  arclength  s.  The  fundamental  unproved  assumption  on 
which  the  method  is  based  is  that  x'(s)  Is  a  continuous  Markoff  process. 

This  means  that  the  probability  that  the  ray  has  a  given  direction  at  s  +  ds 
depends  only  upon  the  probability  distribution  at  s.  Let  us  denote  this  pro¬ 
bability  density  by  P(e,^, s).  Then  as  a  consequence  of  conservation  of  pro¬ 
bability,  P  satisfies  the  Chapman-Kolmogoroff  equation.  The  differential 
form  of  this  equation  is  the  Fokker-Planck  equation 


P  = 
s 


HTe  <6ln  9  Vs  +  -V 

sin  9  r 


lim 


As 


<q2> 

4  As  ’ 


(^5) 


Here  <a  >  is  the  mean  square  value  of  the  angle  a  between  x'(s)  and  x'(s+As) 
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when  x'(s)  is  given.  We  have  calculated  it  in  the  preceding  section  and  it 

2  2 

is  given  by  (4o) .  From  (4o)  it  follows  at  once  that  <sx  >  is  of  order  (As  ) 
for  small  As  and  thus  the  limit  in  (45)  is  zero.  However,  if  we  use  (44)  for 
<x2>  instead  of  (4o),  and  neglect  O(e^),  the  limit  in  (45)  is  D  and  (45)  be¬ 


comes 


P  = 


D 


(sin  9  Pj.  + 


(46) 


sin  0  ^ 

It  would  seem  to  be  inconsistent  to  use  (44),  which  is  valid  for  s 

large  compared  to  the  correlation  distance  a,  and  then  to  let  s  tend  to  zero. 

But  this  kind  of  procedure  is  common  in  applied  mathematics.  It  can  usually 

be  Justified  by  interpreting  it  as  determining  an  asymptotic  expansion  with 

respect  to  some  appropriate  parameter.  In  the  present  case  this  is  so  and 

the  parameter  is  the  correlation  distance  a.  To  show  this  we  shall  compute 

the  limit  of  s  a<xx  >  as  a  tends  to  zero,  by  means  of  (4o).  In  doing  so  we 

-1 

first  observe  that,  from  dimensional  considerations,  U  is  a  function  of  a  r, 
say  N(r)  =  M(a_1r).  Now  (4o)  yields 


—1  2  2  2  -I 

lim  s  a  <a  >  =  lim  4c  <jj.  >  s  a 


a-+o 


a-K3 


°  i 

M(a  J's)-l-s  J'  r  ^"a  Sd'(a  ^r)dr 


=  4* 


<4l2>  s'1 


lim 
a -+0 


-1 
a  s 


aM(a  a)-a-s 


J  x-1M'(x) 


dx 


-  4e2  <n2>  J  x_1M'(x)dx  .  (47) 

o 

2 

From  (47)  and  (42)  we  see  that  for  any  s  ^  0,  (4s)  <a  >  ~  D  as  a  tends  to 

Then  the  limit  in  (45)  is  D.  We  now  see  that  (46)  is  valid  in  the 


zero. 
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limit  in  which  the  correlation  distance  a  tends  to  zero.  Thus  the  solution 
of  (46)  is  the  leading  term  in  the  expansion  of  P  for  small  values  of  a. 

We  have  seen  that  the  assumption  that  x' (s)  is  a  continuous  Markoff  pro¬ 
cess  leads  to  (45).  When  a  tends  to  zero  (45)  becomes  (46).  If  P  is  inde¬ 
pendent  of  which  is  the  case  for  all  s  if  it  is  so  for  s  =  0,  (46)  further 

simplifies  to 


P 


_  D 
s  “  sin  9 


(sin  e  Pe)e  . 


(It8) 


Let  us  now  solve  (48)  subject  to  the  'initial'  condition 

P(0,G,f*)  -  6(9)  •  (49) 

This  'initial'  condition  corresponds  to  the  ray  at  s  «  0  certainly  pointing 
in  the  direction  9=0.  The  solution  of  (48)  and  (49)  is  expressible  In  terms 
of  Legendre  polynomials  Pn( cos  9)  as 

00 

P(s,e,jlO  =  jj-j  2  (2n+l)Pn(cos  9)e"n(n+l)Ds  .  (50) 

o 


From  (50)  ve  see  that  P(»,9,j^)  &  l/4jt  so  that  for  large  s  all  directions  be¬ 
come  equally  probable.  From  (50)  we  easily  find 
-2Ds 

<cos  9>  =  e  .  (51) 

Let  us  compare  (51)  with  the  result  for  <cos  a>  which  can  be  obtained 
from  (4o).  We  observe  that  a  and  9  are  defined  in  the  same  way  so  that  the 
two  results  should  agree.  However,  they  do  not  agree.  Since  (4o)  was  ob¬ 
tained  by  an  honest  method  and  (51)  by  a  dishonest  one,  we  might  be  tempted  to 
conclude  that  the  dishonest  method  has  yielded  an  incorrect  result.  But  the 
preceding  discussion  in  this  section  shows  that  this  is  not  the  proper  conclu¬ 
sion.  Instead  we  must  recognize  that  the  two  results  have  different  domains 
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of  validity.  The  result  (4o)  is  valid  for  small  values  of  s  while  (51)  is 
valid  for  large  values  of  s.  If  we  expand  (4-1)  for  s  large  and  (51)  for  s 
small  we  obtain  agreement  to  two  terms: 

<cos  a>  2;  1  -  2Ds  ~  <cos  ©>  .  (52) 

This  agreement  shows  that  the  two  domains  of  validity  overlap  slightly.  Of 
course  the  question  of  whether  the  result  (51)  is  valid  at  all,  since  it  is 
based  on  the  Markoffian  assumption,  is  not  answered  by  these  considerations. 

The  result  (51)  can  be  uaed  to  determine  the  mean  square  distance  be¬ 
tween  the  endpoints  of  a  ray  of  length  s.  To  use  it  for  this  purpose  we 
write 

s  s 

<|x(s)|2>  =  <  J  x'(t1)dt1*  J  x'  (tg)dtg  > 

o  o 


s  s 


=  JJ  <k-(t1).x-(t2)>dt1dts 


o  o 


s  s 


=  J J  <cos  9(t1,tg)>  dt-^dt^ 


o  o 


s  s 


s  1 
D"  2D2 


For  small  s  this  becomes 
1 2 


=  JJ  exp(-2D|t1-t2| )  dt2dt£ 

o  0 

G  -  *-2D0  • 


(55) 


< 


-]  tL 

x(s)  >  V 


2  2  3 

s  -  ■=■  Ds^  . 

5 


(54) 


The  first  term  of  this  result  does  agree,  but  the  second  does  not  agree, 

r  ~i  2 

>  for 


vith  (35)>  which  gives  the  previously  calculated  value  of  < 


x(s) 


large  s.  That  result  involves  N^( 0)  in  addition  to  D,  while  all  the  re¬ 
sults  of  the  present  method  involve  only  D.  We  must  conclude  that  the 
domains  of  validity  of  (53)  and  (5*0  do  not  overlap. 

Let  us  finally  compute  <^x(s)  *u^>  by  the  present  method.  This  is 
the  mean  square  value  of  the  distance  which  the  ray  travels  along  its  ini¬ 
tial  direction.  It  is  given  by 

s  s 


<(x(8).5)2>.  Jf  <cos  0(t^)cos  ©(tg)>  dt^dtg  . 


(55) 


O  0 


By  using  (50)  we  finally  obtain 


(56) 


Tne  second  form  of  (56)  holds  for  small  Ds.  itB  first  term  agrees  with 
(37) ,  which  gives  the  same  quantity  computed  by  the  previous  method.  How¬ 
ever,  the  second  term  does  not  agree. 

By  substraction  of  (56)  from  (53)  we  obtain  <p  >,  the  mean  square  trans¬ 
verse  displacement  of  the  endpoint  of  the  ray  from  its  initial  tangent  line. 
It  is  given  by 


<V 


(57) 


The  second  form  of  (57)^  which  holds  for  small  Ds,  agrees  vith  (43),  obtained 
by  the  previous  method. 
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These  calculations  suffice  to  illustrate  this  method  and  to  show  to 
vhat  extent  its  results  agree,  or  disagree,  with  those  of  the  preceding 
honest  method. 


6.  Phase  fluctuations 

In  geometrical  optics  rays  do  not  occur  individually  but  only  in  normal 
congruences.  A  normal  congruence  is  a  two  parameter  family  of  rays,  all  of 
which  are  normal  to  some  smooth  surfai- ' .  Every  surface  normal  to  the  rays  of 
such  a  congruence  is  called  a  wavefront.  With  every  normal  congruence  of 
rays  we  associate  a  complex  valued  field  v(x)  which  we  write  as 

v(x)  =  A(x)  elk^X^.  (58) 


The  possibly  complex  factor  A(x)  is  called  the  amplitude  of  the  field  and 
^(x)  is  called  its  phase.  These  quantities  are  uniquely  defined  in  terms  of 
v  by  the  requirements  that  be  independent  of  the  propagation  constant  k 
and  that  A  depend  upon  k  at  most  through  a  factor  of  a  power  of  k.  The  pro¬ 
pagation  constant  k  is  equal  to  2jt/l  where  X  is  the  wavelength  of  the  field 
in  a  region  of  refractive  index  equal  to  unity.  The  amplitude  A(x)  may  be 
either  a  scalar  or  a  vector.  The  real  time  harmonic  field  V(x,t)  determined 
by  v(x)  is  given  by  V(x,t)  =  Rejv(x)e  ivtJ  .  Here  t  denotes  the  time  and  v 
the  angular  frequency  of  the  field. 

Let  x  =  z(a,  p)  be  the  equation  of  a  wavefront  SQ,  i.e.  a  surface  ortho¬ 
gonal  to  a  normal  congruence  of  rays.  We  shall  suppose  that  the  parameters 
a  and  p  are  so  chosen  that  at  aQ, Pq  the  following  equations  hold 


1,  z 


a 


y 


u  * 
a 


z  > 
a 


u_ 


-1 


V 


zaXzp 


u. 


(59) 
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Here  u  denotes  the  unit  normal  to  So  while  p  and  p£  are  the  principle  radii 
of  SQ  at  ao,PQ.  These  equations  imply  that  at  a  , the  parameters  a  and  p 
measure  arclength  along  the  lines  of  principle  curvature  on  S  . 


We  denote  by  x(a, 0,s)  the  ray  normal  to  SQ  at  a, p. 

Thus  measuring  s 

from  S  ,  we  have 
o 

x(a,p,0)  =  z(a,p) 

(60) 

x'(a,p,0)  =  u(a,p) 

(61) 

The  ray  through  a  given  point  y  is  determined  ty  the  values  of  a,  0  and  s  such 
that 

y  =  x(a, p^s)  .  (62) 


Once  these  values  are  determined,  ve  define  the  phase  /(y)  by 
s  s 

1>G)  =  J  njSc(a,p,t)Jdt  -  s  +  c  J  n|x(a,0,t)Jdt  (63) 

o  o 

We  have  taken  i  to  be  zero  on  S  . 

o 

To  evaluate  <j>  we  must  first  solve  (62)  which  we  shall  do  by  expansion  in 
powers  of  e.  By  using  (15)  we  can  write  (62)  as 


y  =  xQ(a,p,s)  +  e  x£(a,p,s)  +  0(e  ). 
Here  x  (a, 0,s)  is  given  by 


m 


x  »  "z(a,  3)  +  su(a,p) 


(65) 


Let  us  set  e  =*  0  in  (64)  and  denote  the  solution  by  ao>0o,SQ.  Then  by  dif¬ 
ferentiating  (64)  with  respect  to  e  and  setting  c  ■  0  we  obtain 

0  *  Voa(VPo'8o>  +  ^(Wol  +  V;(VV'oi  +  VW’o5  • 


(66) 


By  using  (65)  and  (59)  and  the  fact  that  u.xe  =  0  ve  readily  obtain  from 
(66)  the  results 


a 


-z  • x  (a  ,0  ,s  ) 
a  ev  00  o' 

€  1  +  Vi1  ~ 


(68) 


Pe  = 


-*  -*  .  . 

-z  *x  (a  ,8  ,s  ) 
P  €V  o'^V  O ' 


1  +  Sop2 


1 


(69) 


s6  =  0  . 


(69) 


We  now  use  these  results  in  (63).  Since  p  =  Vp.fx  +  a  x  +  Bx  )  for 

€  '  €  e  oa  t  op' 


6  =  0,  ve  obtain 


s 


My)  =  so  +  €  f  Vn[^0(ao,po,t) 

o 


dt 


/ 


+  e  /  Vp 
o 


x  f  x,  ,  3  . t) 

O  •  0  O' 


z  *x  (s  ) 

j  (t)  -  -A-S-Sj.  (1  +  tp-1)Sa 

1  +  s  p,x  -La 


'oPl 


Z  *X  (s  ) 

0  €  O  / 1  -In-* 

* - n  u  +  tp2  'z 


1  +  SoP2 


dt  +  0(e'>)  . 


(70) 


We  next  use  (ll)  for  x^  in  (70)  and  compute  the  mean  value  of  ^(y) .  A 
lengthy  but  rather  straightforward  calculation  yields 


<^(y)>  -  sq+  e£ 


/  2  8o  ^ 

(71) 


If  and  p2  are  both  finite  then,  for  large  s ,  (71)  becomes 
<^(y)>  s»  s  (l  +  e2<(i2>)  +  0(e5)  . 


(72) 


Thus  the  mean  phase  increases  linearly  with  sq  as  if  the  medium  had  the  con- 

2  2 

stant  index  of  refraction  n*  =  1  +  e  <|i  >. 


The  variance  cr 


</>G) 


is  given  by 


^(y) 


=«  < 


< 


n2 


f^(y)  -  <^(y)> 


S  +  +  ~oT  $  “  <35  +  +  -x-  0  >' 

o  ’c  2  ^ee  o  re  2 


>  +  0(e5) 


2</>+0(e3)  . 


(75) 


From  (70)  and  (73)  we  have 


^(y) 


s  s 
0  0 

"■//  <- 

0  o 


[vy 


xo(t2> 


>  dt1dtg  +  O(e^) 


U 

=  2e2  *4i2>  ^  (sQ-r)N(r)dr  +  o(e^)  . 


W 


For  large  sq  this  becomes 


^(y) 


2 €2<U2>  sQ  /  N(r)dr 


(75) 


In  concluding  our  discussion  of  the  phase  we  should  observe  that  ^  is 
a  solution  of  the  eiconal  equation 

(V?02  =  n2  .  (76) 

It  is  a  solution  which  satisfies  the  boundary  condition 


A  =  0  on  S 
r  o 


(77) 


There  are  two  such  solutions  which  merely  differ  in  sign.  Thus  the  analysis 
of  this  section  has  concerned  the  solution  of  a  boundary  value  problem  for  a 
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random  partial  differential  equation  of  first  order.  The  analysis  was  based 
upon  the  solution  of  the  equation  by  means  of  the  characteristics,  which  are 
the  rays  in  the  present  case. 


7.  Amplitude  fluctuations 

The  amplitude  A(x)  also  satisfies  a  first  order  partial  differential 
equation  in  which  ^(x)  occurs.  For  a  scalar  A  it  is 

V •  ( A^Vj^ )  =  0  .  (78) 

Let  us  seek  the  solution  which  has  the  prescribed  value  A°(a, p)  on  the  sur¬ 
face  Sq  given  by  "x  ="z(a,  p).  Thus 


These  equations  can  be  solved  for  A  by  the  method  of  characteristics.  Because 
of ’the  special  form  of  (78)  the  solution  can  be  obtained  easily  by  integrating 
(78)  over  a  volume.  Now  the  volume  integral  of  the  divergence  of  a  vector  is 
equal  to  the  surface  integral  of  the  normal  component  of  the  vector.  To  sim¬ 
plify  this  surface  integral  «e  choose  as  the  volume  a  region  bounded  by  a  tube 
of  rays,  the  watfsfront  Sq  and  another  wavefront  S.  Since  \/<f>  is  tangent  to  the 
rays,  the  surface  integral  over  the  tube  vanishes.  Thus  only  the  integrals 
over  the  two  portions  of  the  wavefronts  remain  and  we  obtain 


(80) 


In  (80)  S'  and  are  the  two  portions  of  the  wavefronts,  do  and  dOQ  are  the 
elements  of  area,  and  we  have  made  use  of  the  fact  that  is  normal  to  the 
wavefronts,  pointing  inward  on  SQ  and  outward  on  S.  Since  (80)  holds  for  every 
choice  of  S'Q,  and  the  corresponding  S'  determined  by  the  rays,  we  can  conclude 


that 
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dcr 


A  |V^|  =  A  (0)  |^(0)  |  dcr 


(81) 


Here  do-o/do'  is  the  ratio  of  corresponding  area  elements  on  So  and  S,  and  the 
argument  zero  denotes  s  =  0  vhich  corresponds  to  a  point  on  S  .  Since 


|V^|  =  n,  we  obtain  from  (8l)  and  (79)  the  result 

C.  -iV/2  vl/2 

Ag(q,.,3)]  .  A°(q,3)  I  -t^.0-) 


^n(x(a,3,s^y  \ 


\  da  J 


(82) 


To  evaluate  A  we  must  compute  dd/daQ ,  vhich  is  the  Jacobian  of  the  map¬ 
ping  of  Sq  onto  S  by  means  of  rays.  This  mapping  is  given  by  the  function 
x(a, 3,s)  with  s  having  the  constant  value  corresponding  to  S.  If  a  and  0 
are  defined  as  before  at  qlq,  on  SQ,  then  along  the  ray  through  this  point 
we  have 


dcr 

dcr 


|^a(VBo'3)  x  VW'D  • 


(85) 


If  we  use  (6*0  for  x(a,0,s),  (85)  becomes 


dcr 

dcr 


[(l+P^s)^  €X£J  X  [(l+p^s)^  exe^j|  +  0(£2) 


(l+p“1s)(l+Pg1s)n  +  €[(1+Pdxs)(zax  xep)  +  (l+p"xs)(x€aX  Zp)] 


-1 


-l_wd  . .  d 


(l+p“1s)(l+p"1s)  +  e[(l+PdXs)u*(zaX  x£p)  +  (l+p"‘Ls)u*(xeaX'zp)]+0(€£::) 


-1 


-1 


(l+p"1s)(l+p"1s)  +  6[(l+p“1s)z(3*x£(3+  (l+p“1s)za-xeJ+0(€2)  .  (81+) 


Nov  (82)  yields  for  A  the  result 


-  d(z)+u(x)J 


>  +  0(e) . 


(85) 


Here  ve  have  introduced  A  vhlch  is  defined  hy 

o 

AQ(a,p,s)  »  A0(a,p)(l+p"1s)"1//2(l+p"1s)"1/2  . 

To  obtain  A  at  a  fixed  point,  say  y,  ve  must  insert  for  a, 0  and  s  in 

2  2 

(85)  the  solutions  of  (62).  They  are  aQ+  ca^  +  0(a  ),  Pq+  eP£+  0(e  )  and 
2 

so+  o(e  ).  Then  (85)  yields 

A(y)  -  Ao(ao,Bo,80)  ji  -  iju+e'1.)-1  (l+p^s)  \-xccl 


(86) 


-  n(z)+u(x) 


>]  +  ‘^caV  vjj  +0(£2)-  (87) 


From  (87)  ve  find  at  once  that 


<A(y)>-  Ao(%>Po>B0)  +  °^6  ^ 


(88) 


Next  ve  compute  the  variance  of  A(y). 

<A2>  -  <A>2  *  <  Ja(J)  -  aJJ  >  +  0(e5) 


(89) 


We  shall  evaluate  it  in  the  special  case  in  vhich  A  *  A  «=  0.  Then  after  a 

oa  op 


lengthy  calculation,  ve  find  for  s  »  a. 
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2.  2^_  3.2 

o  0  €  <n  X!b  A 

-A  >  -  -A>g  .  -  .x.g  » 


5  V  P-L+8  Pg+a, 


pl+s  p2+s 


+  8 


>1  .  Pl+  P2 

+  - —  -  ■  ■■  ■  ■■ 


(p^s)2  (p1+s)(p2+s)  (p2+s)2 


Here 


(90) 


w  w 

C  -  J  J  (v|)2N(r)dr  «  J  [r“2Nrr-  r"5^ 


dr  . 


(91) 


When 


pl  =■  p2  ®  p,  (90)  simplifies  to 


2  2^_  3.2 

2  2  e  I  s*~  21 

<A  >  -  <A>  =  - - g- —  ^  —  +  sp  +  2p^ 


3(p+s)' 


(92) 


If  the  initial  wavefront  is  a  plane,  p  *  »  and  (92)  becomes 


2  2  2  2  2  3  .2 

<A  >  -  <A>  «  -j  €  <4  >  Csy  A 

3  o 


(93) 


If  the  initial  wavefront  is  a  point,  p  =  0  and  (92)  yields 


Q  P  P  X  O 

<A  >  -  <A>  -  <(i  >  Cs^  A 

15  0 


The  result  (93)  is  the  same  as  that  of  Krasilnikov M  vhile  (94)  is  similar 
to,  but  not  identical  with,  that  of  P.G.  BergmannW. 


(9^) 


8.  The  reduced  wave  equation  in  a  random  medium 

Let  us  consider  the  solution  u(x)  of  the  following  problem 

Au  +  k2n2(x)u  *-8(x) 


lim  jx} 


6u 

8|x| 


iknu  ■  0  . 


(95) 

(96) 
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Equation  (95)  is  the  reduced  wave  equation  with  a  source  term  corresponding 
to  a  source  of  unit  strength  located  at  the  origin  x  -  0.  Equation  (96)  is 
the  radiation  condition  which  asserts  that  the  solution  u  describes  outward¬ 
ly  propagating  waves.  This  problem  is  well  posed  for  a  large  class  of  re- 
fractive  indices  n(x).  Let  us  assume  that  n(x)  has  the  form 


n(x)  -  1  +  €|i(x)  . 


(97) 


Of  course  n,  n  and  therefore  u  also  depend  upon  the  variable  <s  which  ranges 
over  the  probability  space  Cl,  but  we  shall  not  write  it  explicitly. 

It  seems  reasonable  to  suppose  that  u  can  be  represented  as  a  power  se¬ 
ries  in  € 

u( Xj  € )  »  UQ(x)  +  €U1(x)  +  €2U2(x)  +  0(e5)  .  (98) 

Upon  inserting  (97)  and  (98)  into  (95)  and  (96),  and  equating  to  zero  the 
coefficient  of  each  power  of  e,  we  obtain 


Au  +  k  u  s  -  6(x) 
o  o 


lim  x 


x  -*  » 


du 
_ _ c 

3lx| 


iku 


2  2 
AUi  +  k  u1  =  -  2k  ^uq 


(99) 

(100) 

(101) 


lim 


.81x1 


-  iku^-  ik|IUo  )  a  0 


(102) 


2  2  2 

AUg  +  k  Ug  ■  -  k  n  uq  -  21c  ^u^ 


du£ 

lim  |x|  ( - iku„  -  ik(iu  I  «  0 

.8|x|  2  11 

1X1  -  » 


(105) 

(io4) 
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The  solution  of  (99)  and  (100)  is 


u  (x)  - 

°  4«|x| 


(105) 


To  solve  (101)  and  (102)  ve  apply  Green's  theorem  to  the  interior  of  a  large 
sphere  and  assume  that  the  integral  over  the  sphere  tends  to  zero  as  its  ra¬ 


dius  becomes  infinite.  Then  ve  obtain 


k2  r  eikl*-*' 


/  iK.|X-X’ 

- h -  ^(x')u  (x')dx'  . 

X-  x' 


(106) 


In  a  similar  vay  ve  can  obtain  the  solution  u_(x)  of  (105)  and  (104),  vhlch 


u2(x)  - 


k2  r  elkl^'l 


^2(x')u0(x')  +  2u(x' )u1(x' )jdx'  .  (107) 


Let  us  nov  compute  the  mean  value  of  u  vhich  requires  the  calculation  of 
the  three  mean  values  <uq>,  <u1>  and  <ug>.  From  (105),  since  uq  is  indepen¬ 
dent  of  u  and  therefore  of  u,  <uo>  «  uq.  From  (106),  since  u1  is  linear  in  u 
and  since  ve  assume  that  <u>  =  0,  it  follovs  that  <u^>  »  0.  To  compute  <U2> 
from  (107)  ve  require  <*i2>  and  <uu1>.  We  assume  that  <ji2(x)>  *  <[?>  is  a  con¬ 
stant,  and  ve  compute  <nu..>by  using  (106)  as  follows: 


<h(x)u1(x)>  =  J 


ikx-x1 


- -  <h(x)u(x'  )>  u  (x'  )dx' 

X  -  x'  ° 


,2  2. 
k  <u  , 


:  p  ik|x-x'| 

-  /  '  ~  - -  N(|x-x'[)u  (x')dx' 

J  X  -  x' 


.2  2 

k  <u  . 


—  N(r)uQ(x+r)dr  . 


(108) 
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In  (108)  we  have  introduced  the  correlation  coefficient  N(|x-x'|)  - 
•4i(x)n(x'  )>/<|i2>)  which  is  assumed  to  be  a  function  of  the  distance  |x-x'  | 
only. 

The  angular  Integration  in  (108)  on  a  sphere  of  radius  r  can  be  per¬ 
formed  explicitly  by  making  use  of  the  following  easily-proved  mean  value 
theorem  which  is  a  consequence  of  (99) 

*  uq(x)  0  <  r  <  |x| 

ik(r-lxl) 

-  uo(x)  sin  k[x|  - — -  r  >  |x|  .  (109) 


When  (109)  is  inserted  into  (108)  it  becomes 


<h(x)u1(x)>  *  uQ(x)2k<u  > 


J'  eikrN(r)sin  kr  dr 


*—  0 


+  e“iklxlsin  k|x|  J'  e2ikrN(r)dr 


-iu0(x)k<Ji  > 


J  (eaiIr-l)j(r)-  J  (e211l<r-|3l>-1)»(r)dr 


0  i-*| 

lxl 

(no) 

For  |x|  »  a,  which  holds  at  points  many  correlation  lengths  from  the  source, 
(110)  can  be  simplified  to 


<p(x)ui(x)>  a  |3<h  >i0(x)  • 

Here  we  have  introduced  the  constant  3  which  is  defined  by 


(HI) 


P  -  -ik  J  ( e21^1*-!)  N(  r )  dr  . 


(112) 


o-l  - 


Nov  vhen  (ill)  Is  used  in  it,  (107)  yields 


.c.  o  ('  ikpc-x'i 

>(*)>  =  t~-  <\i>(  1+20 )  /  —  — -  uo(x')d3c' 

J  x-x' 


<u, 


(113) 


Upon  collecting  our  results  ve  have 


<u(x)>»  uo(x)  +  €2  <ug(x)>  +  O(e^) 


aik|x-x' 


x-x* 


u(x')dx'  +  O(e^).  (ll4) 


Although  (ll4)  gives  <u>  explicitly,  it  is  advantageous  to  rewrite  it  as  an 

2 

integral  equation  for  <U>.  Since  <u>  differs  from  uq  by  o(e  ),  the  dlffer- 

2  2  4 

ence  between  e  <u>  and  e  uq  is  0(c  )  and  we  may  write  (ll4)  as 


2.  2  2  p  lit  |  x-x  -  |  x 

<u(x)>»  uq(x)  +  (1+23)  J  . . <u(x'  )>  dx'  +  0(e5).  (115) 

If  O(e^)  is  omitted,  this  equation  becomes  an  integral  equation  for  <u(x)>. 

It  is  equivalent  to  the  following  differential  equation  and  radiation  condi¬ 
tion  for  <u> 

no 

(H6) 


ikjx-x' j 


A <u>  +  k2(n*)2<U>  *-B(x) 


lim  |x|  -  ikn*u)  *  0  . 

Wl  J 


(H7) 


x  -+  « 


Here  the  effective  index  n*  is  defined  by 
(n*)2  »  1  +  c2  <U2>( 1  +  23) 


S'  00 

»  1  +  s2  <U2>  -  2ik  J  (e21kr-l)N(r)dr^ 


(118) 
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2 

Our  results  (ll6)  and  (117)  show  that  up  to  and  including  terms  of  order  e  , 
the  mean  value  of  u  propagates  as  if  in  a  medium  of  complex  refractive  index 
n*  given  by  (118).  This  is  the  main  result  of  this  section.  The  imaginary 
part  of  kn*  is  the  attenuation  coefficient  a  for  a  wave  in  this  medium.  From 
(ll8)  it  is  given  by 


a  =  2«2  <ji2>k2j^(l  -  cos  2kr)N(r)dr  . 


(119) 


This  result  agrees  with  that  obtained  previously  by  another  method. 
The  solution  of  (ll6)  and  (117)  is 


ikn*|x| 

<u(x>  a  - - - -  •  (120) 


By  comparing  (120)  with  (105)  we  see  that  <u(x)>  is  given  by  replacing  k  by 
kn*  in  uq(x) .  This  same  rule  applies  to  all  pairs  of  solutions  uq  and  <u> 
of  problems  which  differ  from  (95)/  (96)  or  ( 116 ) ,  (117)  merely  in  the  source 
term.  From  (120)  we  have 


|<u(x)>| 


(121) 


The  method  which  we  have  used  to  derive  (116)  and  (117)  can  also  be  used 

to  determine  <u^(x)>  and  from  it  the  variance  <u2(x)>  -  <u(x)>2.  The  quantity 

2 

<u>  is  often  called  the  coherent  wave  and  <u>  the  coherent  Intensity  while 

the  variance  is  called  the  incoherent  intensity.  Since  it  should  be  clear  how 

2  2 
to  obtain  an  equation  for  <u  >  which  is  correct  through  terms  in  €  ,  we  shall 


not  carry  out  this  calculation. 
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